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Exceptional points associated with non-hermitian opera- 
tors, i.e. operators being non-hermitian for real parameter val- 
ues, are investigated. The specific characteristics of the eigen- 
functions at the exceptional point are worked out. Within 
the domain of real parameters the exceptional points are the 
points where eigenvalues switch from real to complex values. 
These and other results are exemplified by a classical problem 
leading to exceptional points of a non-hermitian matrix. 

PACS numbers: 03.65.Vf, 02.30.-f, 05.45.Mt 



I. INTRODUCTION 

Exceptional points are branch point singularities of 
the spectrum and eigenfunctions, which occur generically 
when a matrix, or for instance a Hamilton operator, is an- 
alytically continued in a parameter on which it depends. 
The term 'Exceptional Points' has been introduced by 
Kato [1]. When a physical problem is formulated by 
Hq + XHi with A being a strength parameter, the spec- 
trum and eigenfunctions - E n (X) and \ip n (X)) — are in 
general analytic functions of A. At certain points in the 
complex A-plane two energy levels coalesce. Such coales- 
cence is not to be confused with a genuine degeneracy, 
since the eigenspace of the two coalescing levels is not 
two but only one dimensional; in fact the corresponding 
eigenvectors also coalesce and there is no two dimensional 
subspace as is the case for a genuine degeneracy. 

If both operators, Hq and Hi, are hermitian, these sin- 
gularities - the exceptional points (EP) - can occur only 
for complex parameter values A. As a consequence, at an 
EP the full problem Ho + XHi is no longer hermitian as 
such, but when dealing with matrices, it is still complex 
symmetric. These cases have been studied in some detail 
[2,3] and here we quote the major results. 

EPs are always found in the vicinity of a level repul- 
sion. Suppose that two levels show avoided level crossing 
when A is varied along the real axis; then the analytic 
continuation into the complex A-plane yields a complex 
conjugate pair of EPs where the two coalescing levels are 
analytically connected by a square root branch point [4] . 
The occurrence of EPs is not restricted to repulsions of 
bound states, a recent paper deals with the repulsion of 
resonant states [5]. Being singularities in the interac- 
tion strength EPs determine the convergence radius of 
approximation schemes in the theory of effective inter- 
actions [6]. Quantum mechanical phase transitions are 
characterised by a multitude or accumulation points of 
EPs [7,8]. 



There are a number of phenomena, where the physi- 
cal effect of an EP has been at least indirectly observed. 
Laser induced ionisation states of atoms [9] are a clear 
manifestation of an EP even though in [9] , it has not been 
analysed as such. A recent theoretical paper [10] shows 
that, for a suitable choice of parameters associated with 
an EP, the only acoustic modes in an absorptive medium 
are circular polarized waves with one specific orientation 
for a given EP. Similarly in optics, experimental obser- 
vations in absorptive media [11] reveal the existence of 
handedness since the stable mode of light propagation is 
either a left or a right circular polarized wave for appro- 
priately chosen parameters. This has been interpreted in 
[12] in terms of EPs. A particular resonant behavior of 
atom waves in crystals of light [13] has been interpreted 
[14] in terms of EPs. While absorption is essential in all 
cases, some situations clearly point to a chiral behavior 
of an EP. In fact, the wave function at the EP has been 
shown to have definitive chiral character [15] and this has 
been experimentally confirmed recently [16]. 

In a previous experiment [17], EPs have been inves- 
tigated in a flat microwave cavity. Major findings have 
been the confirmation of a forth order branch point of 
the coalescing wave functions and - depending on the 
path in the complex A-plane - level avoidance associated 
with width crossing or level crossing with width avoid- 
ance. These results are the consequence of the topological 
structure of Riemann sheets at a branch point [4] . The 
experiment thus showed that this topology is a physical 
reality. 

In the following we carry the analysis further in that 
we investigate the EPs of Hq + XHi when Hq or Hi or 
both are no longer hermitian. As a consequence, even 
for real values of A, the problem Hq + XHi is no longer 
hermitian. This lack of hermiticity is different in nature 
from that discussed above where dissipation - like for in- 
stance in the optical model in nuclear physics [18] - makes 
Hq + XHi non-hermitian for complex A. There is a great 
variety of problems in the literature where either the per- 
turbation or the full Hamiltonian is non-hermitian. Bo- 
son mapping [19], effective interactions [20] and the the 
random phase approximation (RPA) in many body the- 
ory [21] yield non-hermitian operators. More recently 
a wider class of non-hermitian Hamiltonians has been 
proposed to address specific symmetries [22] or transi- 
tional points in specific delocalisation models [23] . These 
suggestions have led to a further thorough study [24] of 
non-hermitian operators. 

The present study is motivated by the classical prob- 
lem of two coupled damped oscillators. It gives rise to 
non-hermitian matrices in a natural way. The problem 
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is stated in the following section. The ensuing general 
treatment of section three yields new insights and spe- 
cial features regarding level repulsion. It is shown that 
the change from a complex to a pure real spectrum of a 
(real) non-hcrmitian matrix under variation of the (real) 
parameter A is due to the occurrence of a real EP. As ex- 
pected, the coalescence at the EP of two complex eigen- 
values into one real eigenvalue (which then bifurcates in 
two real eigenvalues) yields only one eigcnfunction in con- 
trast to the usual two for a genuine degeneracy. A typical 
example is the instability point of the RPA. In addition, 
the pattern of level repulsion is distinctly different from 
that of a hermitian problem: the levels approach each 
other in the form of a cusp and not in a smooth way as 
is the case for hermitian Ho and Hi . These general find- 
ings are illustrated in section four where the example of 
section two is resumed. 

We stress that the present paper focusses upon EPs 
and not on the study of non-hermitian operators as such. 
A summary and suggestion is given in the last section. 



II. TWO COUPLED DAMPED OSCILLATORS 

As a first illustration we consider a simple classical 
case of two damped coupled oscillators in one dimension. 
Denoting by Pi,P2,Qi,Q2 the momenta and spatial coor- 
dinates of two point particles of equal mass the equations 
of motion read for the driven system 




= M 




(1) 



dct \iu — Ml = 



and EPs occur for the complex values lu where 



duj 



det \iu - Ml = 



(4) 



(5) 



is fulfilled simultaneously together with Eq.(4). We 
choose the parameter / as the second variable needed 
to enforce the simultaneous solution of Eqs.(4) and (5) 
and keep the other parameters of M fixed, but of course 
any other preference - like choosing g - would be just 
as good and not alter the essential results. Thus we en- 
counter the problem of finding the EPs of the matrix 
problem 
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Note that Mo and Mi are not symmetric. Before we 
turn to explicit solutions and characteristics of Eqs.(4) 
and (5) we first address the general problem of EPs of 
non-hermitian matrices. 
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where uij — ikj, j = 1,2 are essentially the damped fre- 
quencies without coupling and / and g are the coupling 
spring constant and damping of the coupling, respec- 
tively. The driving force is assumed to be oscillatory with 
one single frequency and acting on each particle with am- 
plitude Cj. Here we are interested only in the stationary 
solution being the solution of the inhomogeneous equa- 
tion which reads 




= (ioj - My 1 



exp(iwi). 



(3) 



Resonances occur for the real values to of the complex 
solutions of the secular equation 



III. GENERAL NON-HERMITIAN CASE 

Like for the hermitian operators the behavior around 
an EP can be described locally by a 2 x 2 matrix. The 
reduction of an N dimensional to a two dimensional prob- 
lem is given below. We always consider a situation where 
the unperturbed problem, denoted by Ho, or ho for the 
two dimensional case, is assumed to be diagonal. At first 
we discuss the two dimensional case and assume that the 
Jordan decomposition of the non-hermitian perturbation 
hi, i.e. hi = SJS^ 1 , yields a diagonal matrix J. We thus 
consider 
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(10) 



Note that for hi to be symmetric we would have 4>i = 02, 
i.e. S would be orthogonal. For convenience we have 
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exploited the freedom to use normalized column vectors 
in S. The two eigenvalues of h are given by 

£i, 2 (A) = i(ei+e 2 + A(u;i+W2)±£>) (11) 
with the discriminant 

D = [{ei - e 2 f + A 2 (a>i - uj 2 f + 

\ 1/2 

2A(ei - e 2 )(wi - w 2 )cos(0i + 2 )sec(0i - 2 ) J ■ (12) 
The two levels coalesce when D = that is for 
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FIG. 1. Real spectrum in arbitrary units as a func- 
tion of A. In terms of Eq.(9) the parameters are 
ei = -l,e 2 = l,wi = -0.2, w 2 = -O.6,0i = -2°, 02 = 45°. 
The spectrum is complex between the EPs at A^ p = 7.3 .. . 
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Note that even when all parameters are real the two 
EPs can now occur on the real axis. It happens when the 
signs of 4>\ and (f> 2 arc different. The implication is that 
the spectrum is no longer real when A lies between X^ p 
and A^p. In Fig.l we display a typical case of a spectrum 
of that nature 1 . We recall that at the EPs, where the real 
spectrum ends or begins, only one eigenfunction exists of 
the two by two matrix problem; its precise form is given 
below. Here we stress a general property of a matrix at 
an EP: the Jordan decomposition of 1i(\ep) = TJT^ 1 
yields a non-diagonal matrix J given by the standard 
form 



J, 



EP 



f E(X EP ) 1 
V E(Xep) 



(14) 



1 in Fig.l of [22] the spectra Ek(N) exhibit manifestations of 
real EPs in the variable N. 



whereas for all points A ^ Xep the matrix J\ is diagonal 
and reads 



Jx =T 



1 h(X)T 

_(E 1 (X) 

E 2 (X) 



(15) 



We recall that T is orthogonal (or unitary) only if h(X) 
is hermitian. 

While it is well known that a non-hcrmitian operator 
can have a non-real spectrum, the deviation from the 
hermitian case has, for the real part of the spectrum, 
distinct consequences for the shape of level repulsions; 
this is exemplified in the following section where typical 
results of the two oscillators introduced in the previous 
section are presented. 

Yet the local behavior at the EP is basically the same 
as for the hermitian case. It is clear from Eqs.(ll) and 
(12) that the two eigenvalues are connected at the square 
root branch points situated at A = A PP just as in the 
hermitian case. The difference arises in the eigenfunction 
at the EP. Recall that the coalescence of two eigenvalues 
at the EP is not to be confused with a true degeneracy in 
that there is only one eigenfunction at the EP. At A PP 
this single and unique eigenfunction (up to a possible 
common factor) turns now out to be 



l*ip> = ±M/S£li> + l2>- 



cot ( 



(16) 



We note that, in contrast to the symmetric case, the left 
hand eigenfunction at the EP (or the eigenfunction of the 
adjoint problem) is now different and reads 



Ep\ 



/ tan a>\ . 



tan 



As a result, the relation 



(17) 



(18) 



still prevails just as in the symmetric case. The basis 
vectors \ j), j = 1, 2 refer to the eigenstates of ho. 

We only mention the special case where either (f>i or <p 2 
assume the value or tt/2: in contrast to the hermitian 
case the confluence of the two EPs does not invoke a 
true degeneracy with two independent eigenvectors even 
though it gives rise - for real parameters - to a real level 
crossing; also Eq.(18) is upheld at such points. 

To summarize: for real matrix elements and different 
signs of the angles, the spectrum is complex between the 
two real EPs. Regarding the wave function the quotient 
of the amplitudes of the two coalescing wave functions 
deviates from that of the hermitian case. The genuinely 
complex superposition of the eigenfunction at the EP re- 
mains, however. For real angles <pj (of equal sign), the 
fixed phase difference of ±7r/2 between the basis states 
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occurs just as in the symmetric case, but the ra- 
tio of the modulus of the amplitudes is in general not 
equal to unity. In addition, as the angles may be com- 
plex, not only the ratio of the modulus but also the phase 
difference can be different from the hermitian case. The 
important point is, however, that Eq.(16) describes, up 
to a common factor, the only possible eigcnmodc at the 
EP. 

The reduction locally of an iV-dimensional problem 
to the appropriate effective two-dimensional problem 
around an EP is, mutatis mutandis, achieved along the 
same lines as for hermitian operators [15]. Owing to their 
vanishing norm the two coalescing eigenfunctions domi- 
nate the complete set of all N normalized eigenfunctions 
in the immediate vicinity of an EP. The expansion of the 
TV-dimensional vector 
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in terms of the complete bi-orthogonal set 

X>k(*)><Xfc(A)| = l X^X EP 

with /3^(A) = (Xk(ty\i>Ep)i contains virtually only two 
terms for A — > \ep- In fact we may write 
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up to a common factor; the two non-zero positions 
are given by the values k, k + 1 for which |xfc(A)) and 
|Xfc+i(A)) coalesce. From Eqs. (19,20) the effective two di- 
mensions for any \ip(X)} becomes obvious within a small 
neighbourhood of Xep- 

We do not discuss cases where the Jordan decomposi- 
tions of ho or hi or both do not yield diagonal but block 
matrices as this does not affect the local behavior at an 
EP. This should not be confused with the fact that in all 
cases an EP of the full problem H + XHi (or ho + Xh\) is 
characterized by a non-diagonal matrix J (see Eq.(14)) 
of its Jordan decomposition. 



IV. EXAMPLES 

While there are various physical reasons to consider 
non-hermitian operators, we here focus on the simple 
mechanical model introduced in section two. Note that 
the model can be easily translated into a corresponding 
electronic setting using two coupled R-L-C circuits. EPs 
can always be found for some complex values of the pair 
(u), /), but a complex value of the spring constant / docs 
not appear physical. This is in contrast to quantum me- 
chanical cases discussed previously [4] where dissipation 



is often described by an effective complex interaction. In 
the classical model we therefore introduce the damping 
term of the coupling denoting its strength by the real 
constant g. For given values of cjj and kj we determine g 
such that an EP occurs at a real value of /. The associ- 
ated two coalescing energies are then complex describing 
a damped oscillation being sustained by the driving force. 

Rc oj 
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FIG. 2. Real and imaginary parts of two repelling levels as 
a function of the coupling damping g. The parameters are 
ui = ll>2 = 10, fci = 0.2, ki = 0.1. The spring constant 
/ = 1.005 is chosen such that an EP occurs at g = 0.00075 
close to the real axis. The units are arbitrary. The inset illus- 
trates schematically a typical level repulsion for a symmetric 
matrix 

The particular model reduces the resultant of Eqs. (4) 
and (5) to a polynomial of fifth order in / which is read- 
ily solved. The symmetry of the model implies that an 
EP at the pair (ujep, Iep) is always associated with an 
EP at (—u)* EP , —/bp); we focus our attention on positive 
/, i.e. a repulsive spring, and the physical requirement 
Ssluep < implying proper damping. Obviously, EPs 
can occur only if either oj\ ^ u>2 or k\ ^ ki or both as 
otherwise a genuine degeneracy is found for / = g = 0. 

To get a good understanding for the EPs we first turn 
our attention to the behavior of the eigenvalues of iM 
in Eq.(2) as functions of / and g; the eigenvalues are 
the solutions of Eq.(4). In Fig. 2 the real and imaginary 
parts of two levels coalescing at an EP very close to the 
real g-axis are plotted versus g using for / a real fixed 
value chosen such that an EP occurs in the vicinity of 
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g = 0.00075. 
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FIG. 3. Real and imaginary parts of two repelling levels as 
a function of the spring constant / for fixed value g — 0.00075. 
The parameters are the same as in Fig. 2 

The level repulsion of the real parts and the expected 
crossing [4] of the imaginary parts are distinctly different 
in shape from the usual appearance for hcrmitian ma- 
trices. The cusp originates from the plain square root 
behavior of the singularity, i.e. the difference between 
the two levels is controlled by <~ \/X — Xep (in Fig. 2 
A = g); this is in contrast to the two complex con- 
jugate EPs occurring in the hermitian case where this 
difference is controlled by ~ ^/(A — A#p)(A — \* EP ) = 

\/(A — SRAbp) 2 + (3Abp) 2 and hence produces a smooth 
approach. For illustration a typical shape of the latter 
is drawn schematically in the inset. The deviation from 
a hermitian case is even more dramatic when the two 
levels are plotted against / for a fixed g = 0.00075 as il- 
lustrated in Fig. 3. Yet, the pattern is understood by the 
same mechanism being a square root branch cut running 
along the real /-axis and having a branch point at / s=s 1 . 

As wc deal with a classical system we now turn to the 
behavior of the complex amplitudes q\{uj) and q 2 {tj) of 
Eq.(3). The overall oscillatory time behavior is of no 
interest, we rather concentrate on the modulus and the 
phase difference. In general these complex amplitudes de- 
pend on the amplitudes Cj of the driving force by Eq.(3). 
However, as discussed in the previous section, at the EP 



there is only one mode possible given by Eq.(16) up to 
a global constant factor. In other words, at the EP the 
ratio of the amplitudes of the two coalescing modes is 
given by iyj cot <j)\/ cot 02 which is a function of only Mo 
and M\ and is independent of a driving force, i.e. of the 
Cj. At close distance to ojep the correct value for the 
ratio must therefore approximately be attained. 
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FIG. 4. Top: Modulus of the amplitudes qi(ui) (solid line) 
and q2 (w) (dotted line) for different driving amplitudes (see 
text) ; \qi (ui) | agrees perfectly with \q? (lo) | on the top left. Bot- 
tom: the respective phases indicated in degrees. The dashed 
line is the phase difference between the amplitudes, for conve- 
nience the negative value is plotted. The parameters are the 
same as in Fig. 2 

This is well demonstrated in Fig. 4 where the moduli 
and phases are plotted against the driving frequency for 
the same parameters as in the previous figures giving rise 
to an EP at uiep = 10.05 — 0.15 i. The top drawings are 
the moduli of the amplitudes with the left one referring to 
Ci = i, c 2 = 1 and the right one to c\ — —i, c 2 — 1; below 
are the respective phases. The former choice (left column 
in Fig. 4) is driving the two masses with equal strength 
but with a leading phase of n/2 for the first mass. From 
the drawing we see that qi{u) / q 2 (uo) w +i through the 
whole resonance. This value is almost equal to the exact 
value being qi{ix>Ep)/q2{^Ep) = 0.0049 + i 1.000 .. . (sec 
(21) below). If, however, the 'incorrect' input is enforced 
like the lagging phase (right column), there is more vari- 
ation in the response. Yet, at the resonance the phase 
difference still is +7r/2, i.e opposite to the driving force 
and in line with the mode at the EP, even though the 
ratio of the moduli is quite different from unity. If the 
driving force is getting closer to ll>ep, i-c. if a slightly 
damped excitation is used, the ratio does approach the 
exact value irrespective of the values Cj. In fact, after 
some slightly tedious but straightforward algebra the re- 
sult 

gl {u E p) _ Pi{uep) _ 
92 (Wfip) P2{UEP) 



Jep + ^-2i(g + k 2 )u)EP - u 



EP 



f E p - ZiguEP 
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(ep - Zig^EP , s 

f E p + u\ - 2i(g + h)ujEP - lu 2 ep 

is obtained yielding the numerical value indicated above 
for the parameters considered. While this result is ob- 
tainable analytically for the particular case of Eq.(6), in 
general one has to resort to the two dimensional reduc- 
tion by numerical means and then use Eqs.(9), (10) and 
(16) to find the amplitude ratio. 

V. SUMMARY 

The physical relevance of EPs and their observability 
has been presented in the introduction. Consideration 
of a simple mechanical problem leads to non-hermitian 
matrices and the study of the associated EPs has pro- 
duced new general insights. The parameters for which 
a real spectrum switches to complex values is clearly re- 
lated to the occurrence of EPs on the real axis. The 
instability point of the RPA is just one case in point. 
In addition, the specific shape of the spectrum can be 
quite different from the one encountered for symmetric 
matrices. On the other hand, the topological structure, 
i.e. the Riemann sheet structure of the energy surfaces 
is independent of whether H and/or Hi are hermitian 
or not. The cigenfunctions at the EP have a structure 
similar to the symmetric case except for the value of the 
ratio of the two relevant stat es. This chan ges from ±i 
for the symmetric case to ±ii/cot0i/ cot 2 for the non- 
hermitian case. Note also, that this ratio is different for 
the left hand eigenfunction at the EP. 

The universal significance of the EPs is once more un- 
derlined by the particular example from classical physics. 
While the general features of EPs for non-hermitian Hq 
and H\ have been presented in section three, we be- 
lieve that the particular results of section four can be 
experimentally confirmed, results whose analogues have 
so far been implemented in sophisticated microwave cav- 
ities [16,17], optical systems [11] and atomic spectra [9]. 
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